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The spectral properties of the spinless fermion model with nearest-neighbor repulsive in- 
teractions on a one-dimensional lattice are investigated using the Bethe ansatz. Although its 
bulk quantities are exactly the same as those of the spin- 1/2 XXZ chain, the difference in the 
statistics of particles causes substantial effects on spectral features, such as gapless points of 
dispersion relations and line shapes of spectral functions. In this Letter, we clarify the origin 
of the differences in spectral features between fermionic and bosonic systems in terms of Bethe 
ansatz solutions. We also confirm that the 2-string solutions have considerable spectral weights 
in the high-energy regime. 
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Introduction - In quantum many-body systems, statis- 
tics of particles plays a crucial role in understanding spec- 
tral properties. In conventional fermionic systems in di- 
mensions higher than one, it is known that excitations 
are well understood in terms of fermionic quasiparticles 
in the Fermi liquid theory, where single-particle excita- 
tions are created from a dense distribution of the quasi- 
particles in the region enclosed by the well-defined Fermi 
surface, and the excitations are gapless on the Fermi sur- 
face. On the other hand, in conventional bosonic systems, 
elementary excitations are generated from Bosc-Einstein 
condensates or ordered ground states, where the gapless 
points are determined by the ordering momentum. 

In one-dimensional (ID) chains, the statistics of parti- 
cles may seem to be less important than that in higher di- 
mensions. Indeed, ID hard-core boson models (or equiv- 
alently S=l/2 spin chains) can be mapped onto ID spin- 
less fermion (SF) models through the Jordan- Wigner 
(JW) transformation, 1 ' and bulk quantities are exactly 
the same in these models. However, as will be shown in 
this Letter, the statistics causes substantial differences 
in spectral features, such as dispersion relations and line 
shapes of spectral functions, even for ID chains in the en- 
tire energy range, which naturally connect to low-energy 
properties studied in detail using the field theory. 2 ~ 5 ' 

In addition, systems with repulsive interactions exhibit 
interesting spectral features in the high-energy regime. 
For the spin-1/2 Heisenberg chain, it has been found that 
the high-energy continuum originating from 2-string so- 
lutions of the Bethe ansatz 6 ' has a substantial spectral 
weight in a magnetic field. 7 ' Since the Bethe equation 
of the SF model can be formally written in the same 
form as that of the XXZ model, 1 ' 6 ' we expect that the 2- 
string solutions will also contribute significantly to spec- 
tral weights in the high-energy regime for the SF model. 
Besides, numerical calculations using the time-dependent 
density-matrix renormalization group (tDMRG) method 
have indicated the continuum in the energy and momen- 
tum regions of the 2-string solutions for the SF model. 5 ' 



In this Letter, we confirm the high-energy continuum 
with a considerable spectral weight in the SF model by 
directly using the 2-string solutions. 

Model - In this Letter, we consider the ID SF model 
with the nearest-neighbor interaction: 
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where Cj is the annihilation operator of a fermion at site i, 
and n, is the number operator defined as rij = c\ci. Here, 
t, V, /i, and L represent the transfer integral, interaction 
strength, chemical potential, and number of sites, respec- 
tively. We impose the periodic [antiperiodic] boundary 
condition for chains with an odd [even] number of parti- 
cles in the ground state. 

The ID SF model can be mapped onto the ID spin-1/2 
XXZ model under the periodic boundary condition 
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In this transformation, the fermions in the SF model are 
mapped onto the down spins in the XXZ model with 
J xy = 2t, J z = V, and H = fj,. Hereafter, the number 
of fermions in the SF model (or equivalently the number 
of down spins in the XXZ model) is denoted by M, and 
we assume that L is even. In this Letter, we focus our 
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attention on the case of V/t = 2, which corresponds to 
the Heisenberg chain in a magnetic field (J—J X y=Jz)- 

Method - The eigenstates of the ID SF model can be 
obtained using the Bethe ansatz. 6 ' The wavefunctions for 
V/t = 2 are expressed in terms of rapidities Aj (j = 1 ~ 
M) that satisfy the following Bethe equation: 
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where Ij (j = 1 ~ M) are called Bethe quantum num- 
bers. The difference in the statistics of particles is re- 
flected in these quantum numbers. For the SF model, 
they are integers [half-odd integers] under the periodic 
[antiperiodic] boundary condition. On the other hand, 
for the Heisenberg model, they are integers [half-odd in- 
tegers] under the periodic boundary condition in chains 
with an odd [even] number of down spins. Once a set 
of {Ij} is given, a set of {Aj} is obtained through 
cq. (1). Thus, eigenstates are characterized by distri- 
butions of {Ij}. Distributions of {Ij} are somewhat 
analogous to momentum distributions of noninteracting 
fermions. Namely, {Ij} in the ground state are consecu- 
tively distributed around zero like a Fermi sea, as shown 
in Figs, l(a-f) and l(b-l). Excitations from the ground 
state can be obtained by creating holes [particles] in- 
side [outside] the consecutive distribution, as indicated 
by blue open squares [red solid diamonds] in Fig. 1. The 
hole and the particle are called psinon and antipsinon, 
and denoted by ip and ip*, respectively. 8 ' 

The momentum K and the energy E can be expressed 
in terms of {Ij} or {Aj} as 
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The spectral functions A ± (k,u>) defined below 9 ) can also 
be expressed in terms of {Aj}. 10 ' 



A + (k,oj) = J2\(k^H\GS)\ S(u-€i), 

i 
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where |GS) and \k,€i) denote the ground state and the 
excited state with the excitation energy e; and mo- 
mentum k measured from those of the ground state. 
The spectral functions can be probed by angle-resolved 
photoemission experiments. The A ± {k,uj) correspond 
to dynamical structure factors S ±T (k,u>) for the XXZ 
model, 7 ' 8, n ~ 15 ' where ct and c_fc in eq. (4) are replaced 
by and S£ , respectively, and S ^(k,oj) is reversed 
with respect to uj—0 from that in the conventional defini- 
tion with 5(u>— ej). Note that, even though wavefunctions 
have the same form, |(jfc, e^IGS}! 2 ^ |(fc, |GS)| 2 
because of the difference in the statistics. 10 ' Meanwhile, 
\(k,ei\(n k - 1/2)|GS)| 2 = |<fc, e^IGS)! 2 , 7 ' 8 ' 11 " 17 ) be- 
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Fig. 1. (Color online) Typical distributions of Bethe quantum 
numbers {Ij} in L = 24 for (a) the spinless fermion model and 
(b) the Heisenberg model. Solid circles and red solid diamonds 
denote Ij [j = 1 ~ M). Open circles and blue open squares 
denote unoccupied regions of {Ij}. Blue open squares and red 
solid diamonds behave as quasiparticles called psinons (ip's) and 
antipsinons (ip*'s), respectively. 8 ' The top rows show the distri- 
butions of {Ij} in the ground states with M=8. The second rows 
are those of the lowest excited states with M=9 for < k < 7r. 
The third rows are those of dynamically dominant excitations 
with M=9. The fourth rows are those of the lowest excited states 
with M=7 for < k < w. The fifth rows are those of dynamically 
dominant excitations with M=7. 



cause the number operator is bosonic. 

We calculated A ± (k,uj) using the Bethe ansatz solu- 
tions with real Aj (J = 1 ~ M) of up to 2ip2ip* exci- 
tations. For A + (k,uj), we also considered solutions with 
a 2-string, i.e., a pair of complex rapidities having the 
same real part, 6,18 ' with lipltp* excitations for the real 
Aj (j = 1 ~ M — 2). We calculated /x in the thermody- 
namic limit 19 -' for E in eq. (3). 

Comparisons between fermionic and bosonic systems - 
The results for the spectral functions A ± (k,u>) in the 
SF model and those for dynamical structure factors 
S ±T (k,uj) in the Heisenberg model 7 ' are shown in Figs. 
2(a) and 2(b), respectively. One of the differences is the 
position of the gapless points. In the SF model, the gap- 
less points are located at k — 7r ± kp with the Fermi mo- 
mentum kp = ttM/L as denoted by the vertical dashed 
lines in Fig. 2(a), whereas the gapless point at k = tt 
does not shift as a function of M or the magnetization 
(S z = L/2 - M) in the Heisenberg model [Fig. 2(b)]. 
Another difference appears in the line shapes of A ± (k, ui) 
and S' ±=F (fc, cj). Typical behaviors near the gapless points 
are shown in Fig. 2(d). In the SF model, the peak is very 
sharp and looks almost <5-functional, as shown by the red 
solid line in Fig. 2(d), while that in the Heisenberg model 
has a larger tail in a wide range of energies, as shown by 
the blue dotted line in Fig. 2(d). These differences origi- 
nate from distributions of {Ij} that reflect the statistics 
of particles, as shown below. 
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Fig. 2. (Color online) Spectral functions A^(k,ui) of the spinless 
fermion (SF) model for V/t=2 and dynamical structure factors 
S ^ (fc, u) of the Heisenberg model (J=J xy =J z ). Panels in (a), 
(b), and (c) show results in L=320 for M/L=3/8, 1/4, and 1/8 
(from the left), respectively, (a) A^(k,ui) of the SF model, (b) 
S ±= F(A:,cj) of the Heisenberg model. 7 ' Here, S '"(fc, u>) are re- 
versed with respect to uj=0. The horizontal dotted lines in (a) 
and (b) denote the zero energy, and the vertical lines in (a) de- 
note k = it — kp. The data in (a) and (b) are broadened in 
a Lorentzian form with full width at half maximum 0.08J. (c) 
Contributions to the sum rule for the SF model. Green hatched 
regions denote the contributions from the real {Aj } solutions of 
up to 2tp2ip* excitations for A~(k,uj). Red dotted lines indicate 
those for A + (k,u>). Blue dashed-dotted lines show contributions 
of 2-string solutions with Xtplip* excitations for A+(fc,w). Black 
solid lines indicate the sum of such contributions, (d) Line shapes 
of A+ (fc, to) and S~* (fc, u>) near the gapless points for M/L=3/8. 
The red solid line denotes the line shape of 2ip* excitations in 
the SF model at fc ~ tt — kp = 5tt/S, and the blue dotted line 
denotes that in the Heisenberg model at fc ~ n in L=2240. (e) 
A ± (k,u)) of the SF model for M/L=l/4 in L=320. The green 
dotted line denotes the zero energy, and the red solid line shows 
the dispersion relation of the noninteracting system (V=0). 



Gapless points - In the SF model, when {Ij} in the 
ground state are half-odd integers [integers], those in 
excited states with one fermion added or removed are 
also half-odd integers [integers], as shown in Fig. 1(a). 
Hence, the lowest excited state contributing to A + (k,ui) 
has the distribution of {Ij} with one ip* attached to the 
right [left] edge of the consecutive distribution in the 



ground state as in Fig. l(a-2), which has the momentum 
fc=7r— Uf [ir+kp] measured from that of the ground state 
(eq. (2)). Similarly, the lowest excited state (e»— >+0) in 
A~(k,uj) has the distribution of {/■;} with one ip created 
at the left [right] edge of the consecutive distribution in 
the ground state as in Fig. l(a-4), which has k=n—kp 
[tt+Icf}- Thus, the gapless points are located at k=n±kp 
for the SF model as in the noninteracting case (1^ = 0). 

On the other hand, in the bosonic model, {Ij} in ex- 
cited states with one down spin added or removed are in- 
tegers [half-odd integers] , when those in the ground state 
are half-odd integers [integers] as in Fig. 1(b). Hence, the 
lowest excited state (e,— >+0) in (k, co) and that in 
S ^(k,io) have {Ij} symmetrically distributed around 
zero as in Figs. l(b-2) and l(b-4). As a result, they have 
the momentum k=ir measured from that of the ground 
state (eq. (2)) regardless of the magnetization value. 

Thus, the gapless points differ, depending on the 
statistics of particles. The similarity between excitation 
spectra in the fermionic system (Fig. 2(a)) and those in 
the bosonic system (Fig. 2(b)) may be more easily seen 
by reversing A~(k,u>) (uj < 0) in Fig. 2(a) with respect 
to k = tt — kp (the vertical dashed lines in Fig. 2(a)) and 
shifting momenta in Fig. 2(a) by kp to the right. 

Dynamically dominant excitations - More interest- 
ingly, the above difference in {Ij} distribution between 
the fermionic and bosonic systems also causes the dif- 
ference in line shape between A ± (fc,w) and S ±: ^(k,oj). 
In general, excited states with distributions of {Ij} close 
to that of the ground state carry large spectral weights. 
Thus, dynamically dominant excitations for A + (k,ui) in 
the SF model are obtained by shifting ip* in the low- 
est excited state (the red solid diamond in Fig. l(a-2)) 
away from the center, as indicated by the red arrow in 
Fig. l(a-3). These excitations appear in the momentum 
region \k\ < ir — kp, as shown in Fig. 2(a) for u) > 0. Simi- 
larly, dynamically dominant excitations for A~(k,uj) are 
obtained by shifting rp in the lowest excited state (the 
blue open square in Fig. l(a-4)) toward the center, as 
indicated by the blue arrow in Fig. l(a-5). These excita- 
tions appear in the momentum region \k ± tt\ < kp, as 
shown in Fig. 2(a) for ui < 0. 

These behaviors are in contrast to those in the bosonic 
model. Because of the symmetric distribution of {Ij} in 
the lowest excited state for (k,ui) (Fig. l(b-2)), dy- 
namically dominant excitations are obtained by shifting 
2 -0*'s (the red solid diamonds in Fig. l(b-2)) away from 
the center as in Fig. l(b-3). These excitations form the 
low-energy continuum near k = tt in (fc, w) , as shown 
in Fig. 2(b) for ui > O. 7 ) Similarly, because the lowest ex- 
cited state (ej — » +0) in S~ + (k,uj) has the symmetric 
distribution of {Ij} around zero (Fig. l(b-4)), dynami- 
cally dominant excitations are obtained by shifting 2ip's 
(the blue open squares in Fig. l(b-4)) toward the center 
as in Fig. l(b-5). 8 - ) These excitations form the continuum 
near k = tt for uj < 0, as shown in Fig. 2(b). 7 ' 8 ) 

As seen above, dynamically dominant excitations are 
different between the fermionic and bosonic systems ow- 
ing to the difference in {Ij} distribution. In the SF model, 
large spectral weights are carried by the lip* [lip] mode 
with a small tail 20 -' originating from less dominant multi- 
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ipip* excitations, as shown in Figs. 2(a) and 2(d). In con- 
trast, in the bosonic model, spectral weights spread in the 
2ip* [2ip] continuum, which shows a peak with a larger 
tail, as shown in Figs. 2(b) and 2(d). In this sense, ip*'s 
[ip's] may appear more evenly fractionalized in bosonic 
systems than in fermionic systems. 

High-energy continuum - Next, we discuss the high- 
energy continuum in Fig. 2(a). For the Heisenberg model, 
it has been found that the continuum of the 2-string so- 
lutions has a large spectral weight in (/c,cj), 7 ' which 
has resolved a long-standing puzzle on the high-energy 
properties of quasi-lD antiferromagnets in a magnetic 
field. 7 ' 21 ' Because the Bethe equation of the SF model 
can be formally written in the same form as that of the 
Heisenberg model (eq. (1)), it is naturally expected that 
the 2-string solutions in the SF model will also have con- 
siderable spectral weights in the high-energy regime. Be- 
sides, numerical calculations using the tDMRG method 
have indicated the continuum in the energy and momen- 
tum regions of the 2-string solutions for the SF model. 5 ' 
We directly confirmed it by calculating A ± (k,co) using 
the {Aj} of the 2-string solutions for the SF model. Al- 
though the low-energy continuum from real {A^} solu- 
tions remains almost on the cosine dispersion relation of 
the noninteracting system (the red solid line in Fig. 2(e)), 
the high-energy continuum of the 2-string solutions is 
separated from it, as shown in Figs. 2(a) and 2(e). 

We also investigated their contribution to the sum rule, 

/oo 
duj [A + (k,uj) + A-(k,u)] 
-oo 

= ^[\(k,e i \ c {\GS)\ 2 + \(k,e i \c. k \GS)f , 

i 

as shown in Fig. 2(c). The contribution of the 2-string 
solutions is large near k = 0. Indeed, they occupy more 
than 60% of the total spectral weight per momentum 
near k = for M/L=3/8 and 1/4. The missing weight 
would be mainly accounted for by 3-string solutions as 
in the Heisenberg model. 7 ' It should be noted that, in 
the SF model, the momenta of the 2-string solutions are 
shifted by kp(= irM/L) from those of the Heisenberg 
model, because {Ij} for real Aj (j = 1 ~ M — 2) in the 
2-string solutions are integers [half-odd integers] when 
those in the ground state are integers [half-odd integers] 
as in the case of real {A.,} solutions discussed above. 

We expect that such a high-energy continuum with a 
substantial spectral weight will appear in more general 
ID systems with repulsive interactions. In preliminary 
calculations on the SF model with V > and the XXZ 
model with J z > 0, we have confirmed that the high- 
energy continuum of the 2-string solutions has a consid- 
erable spectral weight in both weak-coupling (XY-like) 
and strong-coupling (Ising-like) regimes. The results will 
be presented elsewhere. Also, it should be noted that so- 
lutions with a pair of fermions in the ID Hubbard model 
have energies of 0(U) in the large- U regime, correspond- 
ing to the upper Hubbard band. 22 ' 

Summary - We have investigated the one-particle spec- 
tral functions of the ID spinless fermion model with the 
nearest-neighbor repulsion using Bethe ansatz solutions. 



The main message is that the statistics of particles sub- 
stantially affects dynamical features, such as the gapless 
points and line shapes of spectral functions, even in ID 
chains. In the SF model, the distributions of Bethe quan- 
tum numbers {Ij} in dynamically dominant excitations 
for A + (k,ut) [A~(k,u>)} have one extra ip* [ip], which re- 
sults in the gapless points at k = ir±kp with kp = irM/L 
and sharp peaks in A ± (k,co). In contrast, in the bosonic 
model, because of the symmetric distribution of {Ij} in 
the lowest excited states, 2ip* [2ip] excitations are dy- 
namically dominant. This results in the gapless point at 
k = ir regardless of the magnetization value and gives 
rise to a larger tail in S ±T (k,td) than in A ± (k,w). 

We also confirmed the high-energy continuum with a 
considerable spectral weight in the SF model by directly 
using the 2-string solutions. We expect that such a high- 
energy continuum will be a general feature in ID systems 
with repulsive interactions. The arguments in this Let- 
ter will be applicable to more general ID fermionic and 
bosonic systems with repulsive interactions. 
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